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ABSTRACT 

PSR J1903+0327, a millisecond pulsar in an eccentric (e = 0.44) 95-day orbit with a 
(~ 1 Mq) companion poses a challenge to our understanding of stellar evolution in bi- 
nary and multiple-star systems. Here we describe optical and radio observations which 
rule out most of the scenarios proposed to explain formation of this system. Radio 
timing measurements of three post-Keplerian effects yield the most precise measure- 
ment of the mass of a millisecond pulsar to date: 1.667 ± 0.021 solar masses (99.7% 
confidence limit). This rules out some equations of state for super-dense matter, fur- 
thermore it is consistent with spin-up of the pulsar by mass accretion, as suggested by 
its short spin period and low magnetic field. Optical spectroscopy of a proposed main 
sequence counterpart show that its orbital motion mirrors the pulsar's 95-day orbit; 
being therefore its binary companion. This finding rules out a previously suggested 
scenario which proposes that the system is presently a hierarchical triple. Conven- 
tional binary evolution scenarios predict that, after recycling a neutron star into a 
millisecond pulsar, the binary companion should become a white dwarf and its or- 
bit should be nearly circular. This suggests that if PSR J1903+0327was recycled, its 
present companion was not responsible for it. The optical detection also provides a 
measurement of the systemic radial velocity of the binary; this and the proper motion 
measured from pulsar timing allow the determination of the systemic 3-D velocity 
in the Galaxy. We find that the system is always within 270 pc of the plane of the 
Galaxy, but always more than 3 kpc away from the Galactic centre. Thus an exchange 
interaction in a dense stellar environment (like a globular cluster or the Galactic cen- 
tre) is not likely to be the origin of this system. We suggest that after the supernova 
that formed it, the neutron star was in a tight orbit with a main-sequence star, the 
present companion was a tertiary farther out. The neutron star then accreted mat- 
ter from its evolving inner MS companion, forming a millisecond pulsar. The former 
donor star then disappears, either due to a chaotic 3-body interaction with the outer 
star (caused by the expansion of the inner orbit that necessarily results from mass 
transfer), or in the case of a very compact inner system, due to ablation/accretion 
by the newly formed millisecond pulsar. We discuss in detail the possible evolution of 
such a system before the supernova. 
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1 INTRODUCTION 

PSR J1903+0327 was the first millisecond pulsar (MSFQ) 
discovered in the on going Arecibo L-band Feed Array 
(ALFA) p ulsar survey ICordes et all (120061 ) . In the discov- 
ery paper IChampion et al.l (|2008l ). we presented the results 
of phase-coherent radio timing of this pulsar carried out with 
the Green Bank and Arecibo radio telescopes. These quickly 
revealed that the pulsar was in a 95-day orbit around a 1 
solar mass (Mq) companion. This object is remarkable for 
being the first (and thus far, the only) disk MSP known 
to have an eccentric (e = 0.44) orbit. In globular clusters 
(GCs) there are several binary MSPs with eccentric orbits; 
but those are thought to be caused by perturbations of the 
binary systems by occasional close interactions with other 
stars. 

Coincident with the pulsar position derived from the 
timing, a star was found whose near-infrared magnitudes 
were consistent with a 1 Mq main-sequence star at the dis- 
tance and reddening estimated for PSR J 1903+0327. It was 
not known then whether this was just an unlikely (~ 2.6%) 
chance alignment or whether the star is genuinely associ- 
ated with PSR J1903+0327, and if so whether it is the 
binary companion responsible for the 95-day orbit of the 
pulsar. Such a finding would be surprising, as the conven- 
tional understanding of MSP evolution posits that such a 
neutron star (NS) is spun up to high spin frequencies by 
accretion of matter and angular momentum from a com- 
panion _stax_wdiite_th«_TOmr^ a giant 
phase iBhattacharva fc van den Heuvell l|l99ll ): this circu- 
larises the system, and a recycled MSP is left orbiting a 
low-mass white dwarf (the remnant core of the donor) in 
a low-eccentricity orbit (e < 10~ 3 ; Phinney 1992). Until 
the discovery of PSR J1903+0327 all known MSPs in the 
Galactic disk had such low-eccentricity orbits. For reviews, 
see Phinney & Kulkarni (1994), Stairs (2004), Tauris & van 
den Heuvel (2006). 

For these reasons, Champion et al. (2008) proposed that 
PSR J1903+0327 may be part of a triple system where 
the 95-day orbit of the pulsar is caused by a massive un- 
seen WD and the third member is the star detected in 
the near-infrared. The latter is in a long-period orbit and 
drives the eccentricity o f the inner pair through the Kozai 
mechanis m iKozail Jl962l). An alte rnative possibility, also dis- 
cussed in IChampion et al] l|200ct ) is that the companion to 
PSR J1903+0327 in the 95-day orbit is the star detected 
in the near-infrared, but that this eccentric, unusual sys- 
tem originated in an exchange interaction in a dense stellar 
environment, like a globular cluster. 

In this paper, we present new optical measurements and 
further radio timing of PSR J1903+0327 obtained with the 
aim of testing these scenarios. The plan for the rest of this 
paper is as follows. The optical and radio observations are 
described §[2] The immediate results from these observations 
are described in § [3] In § [4] we discuss the implications of 
these results regarding the formation and evolution of this 
system. In §[5]we discuss how this system might have formed. 
We summarise our main conclusions in § [6] 



1 We define a millisecond pulsar as a pulsar with spin period 
P < 20 ms and with a low surface magnetic field, B ~ 10 8-9 G. 
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Figure 1. This 32" X 32" subsection of a 5min SDSS i-band 
image taken with GMOS at Gemini North on Sept. 20th, 2007 
shows the location of the 1" slit. Besides the companion to 
PSR J1903+0327 denoted with 'PSR', spectra were also ex- 
tracted of stars A, B and C. 



2 OBSERVATIONS 
2.1 Optical observations 

Long-slit spectroscopy of the suspected counterpart 
to PSR J1903 + 0327 was obtained with FORS2 
lAppenzeller et all (|l998T ). the low dispersion spectro- 
graph of ESO's Very Large Telescope. Four spectra were 
obtained on 2008 June 21, three on 2008 August 23 and 
one a day later, on August 24. All spectra had exposure 
times of 46 minutes, and used a 1" slit combined with the 
1028Z holographic grism, providing wavelength coverage 
over 7830 A to 9570 A. The detectors were read out with 
2x2 binning, yielding a resolution of 3.4 A, sampled at 
0.86 A pix -1 . The slit was placed such that both the pulsar 
companion and a bright nearby star were centred on the slit. 
The observations were taken during clear and photometric 
nights, with the seeing between 0'.'48 and 0'/72. The spectral 
observations were corrected for bias and flat-fielded using 
lamp flats. 

Spectral extraction is complicated by the bright star, 
henceforth star A, located 2 / .'3 from the pulsar counterpart 
(see Fig. 12. 1[) . The star is brighter by about 5.6 mag in the 
/-band and as a result about 20% of the detected counts at 
the spatial position of the pulsar counterpart are from the 
wings of the brighter star. The regular optimal extraction 
algorithm as described by Home (1986) is not suited to ex- 
tract blended spectra as it makes no assumptions about the 
spatial profile. Instead, we use a variation on the algorithm 
by Hynes (2002), which is based on the optimal extraction 
method of Home (1986). 

Hynes (2002) uses an analytic function to describe the 
spatial profile as a function of wavelength. This profile is fit- 
ted to the spectrum of an isolated template source and then 
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Figure 2. The spectrum of the counterpart of PSR J1903+0327. Shown with error bars are the average data from the 8 individual 
spectra, normalised, shifted to zero velocity and binned by a factor of 2. The best fit normalised spectrum from the library of Munari et 
al. (2005) is shown as the solid black line, shifted downwards by 0.5 units. The average night sky emission spectrum is shown by the solid 
grey line, showing the location of emission lines as regions of reduced sensitivity. The vertical scale is arbitrary. The location of stellar 
absorption lines belonging to the Ca II triplet and the Paschen series of Hydrogen are denoted in grey above the spectrum. 



used to simultaneously extract the spectra of the blended 
sources. We use a variation on this method. First, instead of 
using a Voigt function to describe the spatial profile we use 
a Moffat (1969) function, essentially a modified Lorentzian 
with a variable exponent. Secondly, instead of removing the 
sky contribution before extraction, we include it in the fit, 
representing the sky as a first order polynomial added to 
the Moffat profiles for each object in the blend. Finally, the 
absence of an isolated template source forced us to use star 
A as a template. We used an iterative scheme to converge 
the properties of the profile as a function of wavelength from 
star A while removing the contribution of the pulsar coun- 
terpart and another faint star (star B) that was also part of 
the blend. In addition to these three objects, a fourth star 
(star C) is located 28" North- West of PSR J1903+0327 on 
the slit. This star was not blended and extracted normally. 

Arc lamp exposures obtained during daytime with the 
telescope pointing towards the zenith were used for wave- 
length calibration. However, comparison of wavelengths of 
the night sky emission lines from the science exposures 
showed wavelength offsets of up to 0.8 A between different 
exposures, most likely due to flexure caused by the telescope 
pointing away from zenith. To correct for these offsets the 
wavelength calibration of June 21st was used to create a sec- 
ondary line list of some 60 night sky emission lines in the 
first spectrum, which was subsequently used to calibrate the 
remaining 7 spectra. Typical rms residuals of these fits were 
less than 0.05 A. 

Figure [2] shows the averaged spectrum of the counter- 
part to PSR J1903+0327. The spectrum has a signal-to- 
noise ratio (S/N) of 15 to 20 in the wavelength range of 
8000 A to 9000 A. The spectrum is generally featureless ex- 
cept for clear absorption lines around 8540 A and 8660 A. 
At these wavelengths the most common spectral lines are 



those belonging to the Call triplet, a nd the Paschen series 
of Hydrogen (see ICenarro et al.ll200ll ). Though the Paschen 
series has lines on these wavelengths, the absence of the 
other lines in the series argues for them being due to the 
Calcium triplet. 

2.2 Radio timing observations 

We made pulse time of arrival (TOA) measurements of 
PSR J1903+0327 using the 305-m Arecibo radio telescope 
and the 105-m Green Bank Telescope from 2006 Decem- 
ber through 2010 January; the observation dates and radio 
frequencies are summarised in Fig. [3] The TOAs obtained 
before 2008 January were described and used by Champion 
et al. (2008). We now describe the observational setups used 
at the two observatories. 

In this re-analysis we use all the TOAs obtained from 
processing of data taken with the Green Bank S-band re- 
ceiver (with frequency coverage from 1650 to 2250 MHz); 
these cover most of the year 2007. The data were acquired 
using the "Spigot" pulsar backend — a 3-level autocorre- 
lation spect rometer which samp les observing bands of up 
to 800 MHz lKaplan et all (|2005l ). Autocorrelation functions 
(ACFs) of length 2048 lags were accumulated with 3-level 
precision and written to disk every 81.92/is. The ACFs 
were subsequently Fourier transformed to synthesise 2048- 
channel power spectra. We separately analyzed the lower 
and upper half of the band, dedispersing each and deriving 
TOAs using the methods described in the Supplement of 
Champion et al. (2008). The TOA uncertainties are 6 /is in 
the upper half and 10 /is in the lower half of the band. In 
Champion et al. (2008), only data from the upper half of 
the band were used. 

The Arecibo data were obtained with four Wideband 
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Figure 3. Top: TOA frequencies versus time (in years). Bot- 
tom: Post-fit TOA residuals versus time (in years). The different 
colours indicate different observing systems. Light Blue: Green 
Bank SPIGOT data, centred at 2100 MHz, Pink: Green Bank 
SPIGOT data, centred at 1800 MHz. This was not previously 
taken into account. Dark blue: Arecibo WAPP data taken with 
the "S-high" receiver. Light Green: Arecibo WAPP data taken 
with the "S-wide" receiver. Red: Arecibo WAPP data taken with 
the "L-wide" receiver. Black: Arecibo ASP data (also taken with 
the "S-wide" receiver). 



Arecibo Pulsar Processors (WAPPs; Dowd, Sisk & Ha- 
gen, 2000). With the "L-wide" receiver, the WAPPs had 
bandwidth 50 MHz and were centred at 1520, 1570, 1620 
and 1670 MHz; each produced 512-lag ACFs sampled with 
3-level precision, accumulated and integrated every 64/is. 
With the "S-wide" receiver, the WAPPs had bandwidth 50 
MHz and were centred at 2125, 2175, 2225, and 2275 MHz, 
though the latter was sometimes disconnected in order to 
use the Arecibo Signal Processor (ASP, described below) in 
its place. At these frequencies, the WAPPs produced 256-lag 
ACFs which were integrated for 32/is. With the 'S-high' re- 
ceiver, the WAPPs were configured in dual-board mode, al- 
lowing 100 MHz of bandwidth in each of 8 correlator boards, 
centred at 100 MHz intervals from 3150 to 3850 MHz. Each 
board produced 128-lag ACFs integrated for 32/is. 



As for the Spigot data, the ACFs from the WAPPs are 
Fourier transformed to produce power spectra. These data 
were then dedispersed and folded modulo the pulsar's pe- 
riod using a routine written by one of us (IHS) specifically 
for this purpose; this folds the spectra using polynomial co- 
efficients calculated specifically for the specified observing 
frequencies. The folding results in 128-bin pulse profiles ev- 
ery 500 seconds for each WAPP. 

During some of the S-wide observations we collected 
data with the Arecibo Signal Processor (ASP, Demorest 
2007) in parallel with three WAPPs. The ASP band had 
a centre frequency of 2350 MHz. The ASP coherently dedis- 
perses a maximum of 16 4-MHz bands, for a maximum band- 
width of 64 MHz. These data were then immediately folded 
at the pulsar's rotational period, and 512-bin pulse pro- 
files were stored for each of these bands every 120 seconds. 
However, because of the large dispersion measure (DM) of 
PSR J1903+0327 we can only coherently dedisperse 12 of 
these 4-MHz bands at any given time. Finally, we add all the 
bands and four 120-s integrations to produce pulse profiles 
with good S/N. 

All puls e profi les are cross-correlated in the Fourier do- 
mam iTavlorl (|l992T ) with a low-noise template derived from 
the sum of the pulse profiles obtained with the same spec- 
trometer and at the same frequency. From this, we derive 
a total of 1872 topocentric TOAs. Both the high-resolution 
WAPP and ASP data provide TOA measurements with an 
average uncertainty of ~ 1/xs. 

In the next step we use the tempo2 software pack- 
age (Hobbs et al. 2006) for TOA analysis. This applies the 
clock corrections intrinsic to the observatory, the Earth ro- 
tation data and the observatory coordinates to convert the 
TOAs to Terrestrial Time (TT), as maintained by the Bu- 
reau International des Poids et Mesures. For conversion of 
TT TOAs to Coordinated Barycentric Time (TCB) TOAs 
we used the DE/LE 421 solar system ephemeris (Folkner, 
Williams & Boggs, 2008). The program then minimises the 
squares of TOA residuals — the difference between the 
observ ed and predicted TOAs. We use d the "DD" orbital 
model iDamour fc Deruelld l| 19851 , ll986T ); this is optimal for 
describing eccentric orbits in a theory-independent manner. 
The best-fit parameters are listed in Table 1, their uncertain- 
ties are the 1-a values estimated by tempo2 (throughout 
this paper, the quoted uncertainties are 1-a, except when 
stated otherwise). The TOA residuals are displayed in the 
bottom plot of Fig. [3] 

We have performed a parallel analysis using the tempo 
software packagfl The resulting timing parameters and 
their uncertainties (as estimated by tempo) are also pre- 
sented in Table 1. Several timing parameters are different 
from those estimated by tempo2. All time-like quantities 
(spin frequency, orbital period, projected semi-major axis, 
time of passage through periastron and uj, which is strongly 
covariant with the latter) are different because tempo con- 
verts the TOAs to Barycentric Dynamic Time (TDB), a 
time scale that is different from TCB. The right ascension 
and declination ar e differen t beca use the DE/LE 405 solar 
system ephemeris IStandishl (|l998T ) uses an earlier celestial 
reference frame. Furthermore, we included in tempo a mod- 
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ified version of the DD orbital model that uses the ortho- 
m etric paramete r izatio n of the Shapiro delay, as described 
in lFreire fc Wexl (|2010l ). The "orthometric amplitude" (h s ) 
and the "orthometric ratio" (?) provide an improved descrip- 
tion of the areas of the (m c ,sini) plane where the system 
is likely to be compared to the "range" (r) and "shape" (s) 
parameters used in the normal parameterization. In cases 
where we can measure other post-Keplerian (PK) parame- 
ters, the new parameterization yields an improved test of 
general relativity. 

Using either tempo or tempo2 there are no apparent 
trends in the residuals and the normalised \ 2 is 1-3, which 
is similar to what we observe in other MSPs timed with 
these observing systems. This means that the timing model 
provides a complete description of the observed TO As, with 
no significant, unmodelled effects present. This also means 
that the \-a uncertainties reported by tempo and tempo2 
are reasonably accurate estimates of the real uncertainties. 
The non-unity value of the normalised \ 2 likely reflects 
a slight underestimation of the TOA uncertainties by the 
cross-correlation analysis, but it could also be due to red 
noise in the pulsar rotation, as suggested by the marginal 
(almost 3-cr) detection of v. To account for this, the timing 
parameters in Table 1 were obtained after the multiplication 
of the TOA uncertainty estimates by a scale factor; this is 
calculated separately for each dataset so that its normalised 
X 2 is 1. This scale factor is 1.1 for all the Arecibo data and 
1.3 for the GBT data. 



3 RESULTS 

3.1 On the nature of the companion 

Radial velocities and spectral parameters were determined 
by comparing the observed spectra, both individual and av- 
erages at each epoch, with synthetic spectral templates from 
the spectral library of Munari et al. (2005). These models 
range in effective temperature T c g, surface gravity g, rota- 
tional velocity V ro t and metallicity [M/H] and are sampled 
at a resolution of A/AA = 20000. The wavelength range be- 
tween 8400 A and 8800 A was used. To remove the effects of 
the finite resolution of the instrument, the synthetic spectra 
were convolved with a Gaussian profile with a width equal 
to the seeing, and truncated at the width of the slit. For 
each object, an iterative procedure was used to first deter- 
mine the radial velocities using a starting template, then use 
those velocities to shift all spectra to zero velocity and create 
an average which was then used to obtain a better matching 
template for determining more accurate radial velocities. 

The comparison of the a verage spectrum with synthetic 
spectra from the library of iMunari et alj (|2005T ) yields a 
temperature of T e g — 5825 ± 200 K (la) and sets a 2-cr 
lower limit on the surface gravity of log g > 4cms~ 2 , firmly 
excluding the possibility of a giant or sub-giant star (see 
Fig. [4]). Based on the shape of the absorption lines we es- 
timate a 3-cr upper limit on the rotational broadening of 
Vrot sin it, < 140 km s -1 , where i* is the inclination of the 
star. The observed spectrum does not constrain th e metallic- 
ity. Stellar evolution models iGirardi et al.l (|2000T l show that 
a main sequence (MS) star with the observed mass and ef- 
fective temperature of the companion of PSR J1903+0327 
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Figure 4. Confidence contours for the effective temperature 
(T e fr) and surface gravity [log(g/cms~ 2 )] for the companion of 
PSR J1903+0327. 
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Figure 5. The predicted radial velocities of the pulsar (solid 
line) and companion (dashed line) based on the pulsar timing 
ephemeris. These are vertically offset by the systemic radial ve- 
locity of the system relative to the solar system barycentre 7, 
which cannot be determined from the radio measurements. The 
two optical radial velocity measurements are shown with their er- 
ror bars. Their difference is consistent with the predictions from 
the pulsar timing ephemeris, and their absolute values allow an 
estimate of 7 = 44.3 ± 4.9 km s _1 . 

is consistent with the observed infrared colours for distances 
between 6 and 8 kpc, assuming the estima ted reddening as 
a func tion of distance from 2MASS stars IChampion et alj 
(2008). This is in agreement with the 6.4 kpc distance 
estimated from the p ulsar's DM (see Table 1) using the 
ICordes fc Laziol i|2002l ) model of the Galactic electron dis- 
tribution. The predicted ages from these models range from 
4 to 7Gyr. 

3.2 Radial Velocities 

Table [5] shows the radial velocities of the pulsar counter- 
part relative to the Solar System barycentre. They are 
Vi = 92.4 ± 8.3 krns" 1 during the first epoch, and V2 = 
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Timing parameters 


Timing Program 


TEMPO 


TEMPO 2 


Solar System Ephcmcris 


DE 405/LE 405 


DE 421/LE 421 


Reference Time Scale 


TDB 


TCB 


Orbital Model 


Modified DD 


DD 


Reference Time (MJD) 


55000 


55000 


Right Ascension, a (J2000) 


19 h 03 m 05!793296(2) 


19 h 03 m 05?793213(10) 


Declination, <5 (J2000) 


03° 27' 19^21053(6) 


03° 27' 19'.'20911(6) 


Proper Motion in a, fi a (mas yr -1 ) 


-2.01(7) 


-2.06(7) 


Proper Motion in S, fig (mas yr -1 ) 


-5.20(12) 


-5.21(12) 


Spin Frequency, v (Hz) 


465.135245551237(10) 


465.135238339217(9) 


First Derivative of v, v (10 Hz s _1 ) 


-4.0719(2) 


-4.0719(2) 


Dispersion Measure, DM (cm~ 3 pc) 


297.5244(6) 


297.5245(6) 


First Derivative of DM (cm - 3 pc yr~ 1 ) 


-0.0083(6) 


-0.0084(6) 


Orbital Period P& (days) 


95.174117277(14) 


95.174118753(14) 


Projected Semi-Major Axis, x (lt-s) 


105.5934628(5) 


105.5934643(5) 


Eccentricity, e 


0.436678410(3) 


0.436678409(3) 


Longitude of Periastron, ui (°) 


141.6531042(2) 


141.6524786(6) 


Time of Passage through Periastron, To (MJD) 


55015.58140451(4) 


55015.58158859(4) 


Derivative of x, x a (10 _15 lt-s s -1 ) 


+20(3) 


+21(3) 


Apsidal Motion, ui ("yr -1 ) 


0.0002400(2) 


0.0002400(2) 


"Range" Parameter of the Shapiro Delay, t/Tq (Mq) .... 




1.03(3) 


"Shape" Parameter of the Shapiro Delay, s 




0.9760(15) 


Orthometric Amplitude of the Shapiro Delay, /13 (/is) 


2.602(25) 




Orthometric Ratio of the Shapiro Delay, f 


0.803(6) 




Limits (not fitted with other timing parameters) 


Second Derivative of v, v (10 26 Hz s 2 ) 


4.8(1.7) 


6.1(1.6) 


First Derivative of e, e (10~ 16 s" 1 ) 


1.4(6) 




First Derivative of P b , P b (lO^ss -1 ) 


-53(33) 


-64(31) 


Derived Parameters 


Spin Period, P (ms) 


2.14991236434921(7) 




First Derivative of Spin Period, P (lO^ss -1 ) 


1.88203(11) 




Characteristic Age, t c (10 9 yr) 


1.9 




Dipolar Magnetic Flux Density at the Poles, Bo (10 8 G) . . . 


2.0 




Galactic Longitude, I 


37? 3363 




Galactic Latitude, b 


-1?0136 




Distance, D (kpc) 


6.4 ±1.0 




Total Proper Motion, fj, (mas yr -1 ) 


5.60(11) 




Galactic Position Angle of Proper Motion, © M 


263.9 ±0.8° 




Velocities, Solar System Baryccntrc Reference Frame 






Transverse velocity, Vt (km s _1 ) 


168 ±24 




Radial velocity, 7 (km s _1 ) 


42.0 ±4.4 




Total 3-D velocity, V (km s" 1 ) 


174 ±25 




Velocities, Galactocentric Reference Frame (cylindrical) 






V p (km s- 1 ) 


17 ±13 




V (km s- 1 ) 


-189 ±5 




V z (km s- 1 ) 


-11±4 




Total 3-D velocity, V (km s" 1 ) 


190 ±5 




Velocities, Relative to Pulsar Standard of Rest (cylindrical) 






V (km s- 1 ) 


28 ±5 




Total 3-D velocity, V (km s" 1 ) 


37±9 




Derived Masses and Orbital Orientation 






Mass Function, / (Mq) 


0.139558441(2) 




Orbital Inclination, i (°) 


77.47(15) or 102.53(15) 




Total Mass of Binary, M t (Mq) 


2.697(29) b 




Companion Mass, m c (Mq) 


1.029(8) 6 




Pulsar Mass, m p (Mq) 


1.667(21) 6 





Table 1. Timing and derived parameters for PSR J1903+0327 using TEMPO and TEMP02 as described in the text. All parameters are 
as measured at the Solar System Barycentre. In parentheses we present the 1-a uncertainties, except in a few cases ( b ) where we present 
99.7% confidence limits. 
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Table 2. The celestial position, JHK B magnitudes and averaged radial velocities (relative to the solar system barycentre) at both epochs 
of the four stars located on the slit. The position and magnitudes are obtained from the data described in Champion et al. (2008). 

J H K s Vi (kms" 1 ) V 2 (kms" 1 ) 

14.85 ±0.09 14.44 ±0.10 14.23 ±0.19 42.0 ± 0.2 56.3 ± 0.2 

19.49 ±0.09 18.52 ±0.10 18.08 ± 0.09 49.6 ±6.1 48.6 ± 2.7 

15.07 ±0.09 14.52 ±0.10 14.27 ± 0.09 45.9 ± 0.2 54.3 ± 0.4 

19.22 ±0.09 18.41 ±0.10 18.03 ± 0.09 92.4 ± 8.3 20.2 ± 5.2 



ID 


Q J2000 


<5j2000 


A 


19 h 03 m 05?751 


+03°27'21'.'31 


B 


19 h 03 m 05?703 


+03°27'23^'83 


C 


19 h 03 m 05?106 


+03°27'47'.'50 


PSR 


19 h 03 m 05?803 


+03° 27' 19'.' 18 



20.2 ± 5.2 kms -1 during the second. The velocity difference 
between the two epochs of V\ — V2 = 72.2 ± 9.8 kms -1 de- 
viates from constant velocity at the 7Aa level. Between the 
two epochs, the averaged radial velocities of stars A and 
C differ by Vi - V 2 = -14.2 ± 0.3kms -1 for star A and 
—8.4 ± 0.4 kms -1 for star C. These differences are caused 
by slight differenc es in centring of the source on the slit 
iBassa et al.l (120061 ). Based on the location of star A with 
respect to the centre of the slit before and after each expo- 
sure, we estimate an average offset in radial velocity of about 
—8.6 ±3.4 kms -1 for the first epoch and 3.7 ± 2.9 km s -1 for 
the second epoch. This is comparable to the weighted aver- 
age in Vi — V2 = —12.1 ± 0.3kms -1 for stars A and C. 

Subtracting these radial velocity offsets from the mea- 
sured velocities of the counterpart, we obtain Vi = 101.0 ± 
8.7kms -1 and Vi = 16.5 ± 5.9 kms -1 , and a velocity dif- 
ference of Vi — V2 = 85 ± 11 kms" 1 . This is consistent with 
the Vi — V2 = 88.72 kms -1 difference in radial velocity pre- 
dicted from the pulsar timing ephemeris. This confirms that 
the optical counterpart to PSR J1903±0327 is the object 
in the 95-day orbit around the pulsar. Fig. [S] shows the ra- 
dial velocity measurements of the companion and the radial 
velocity predictions based on the orbital parameters deter- 
mined from pulsar timing. This provides an independent es- 
timate of the mass ratio (R = 1.55 ±0.20) and the systemic 
radial velocity of the binary, 7 = 44.3 ± 4.9 kms -1 , which 
cannot be derived from the radio timing. 



3.3 Shapiro delay 

Since the companion is a non-degenerate star, it could in 
principle have a strong stellar wind, which should be de- 
tectable as a variation of DM as a function of orbital phase. 
Such a signal would produce a distortion in our measurement 
of the Shapiro delay. In Fig. [6] we display the DM averaged 
over 36 bins of the pulsar's mean anomaly, after correction 
of the long-term DM variations. For some of these intervals 
we have smaller amounts of data or a small range of frequen- 
cies. If these result in DM determinations with uncertainties 
greater than 0.01 cm -3 pc, they are not depicted. We detect 
no DM variations greater than 0.001 cm -3 pc. This means 
that the mass estimates derived from the Shapiro delay are 
accurate, within their uncertainties. 

The better timing precision, larger number of TOAs, 
longer data span and optimised orbital coverage of our new 
dataset not only improve previous measurements of rela- 
tivistic effects (namely the observed apsidal motion Cj and 
the s parameter of the Shapiro delay, see Champion et al. 
2008) but also allow a precise measurement of the r param- 
eter of the Shapiro delay, in addition to a measurement of 
the proper motion /x and the variation it causes on the ap- 



parent size of the orbit x (see Table 1). We now discuss the 
relevance of these parameters. 



3.3.1 Shapiro delay and component masses 

With the previously available timing data it was not pos- 
sible to measure the companion mass m c from the Shapiro 
delay alone. The precise (and frequent) Arecibo timing has 
completely changed this situation, providing a very clear 
Shapiro delay signal (see Fig. 0). Assuming that general rel- 
ativity (GR) is the correct theory of gravity, we obtain 



r 



= 1.03 ±0.03, 



(1) 



where T = GMq / c 3 = 4.925490947/is is the solar mass M© 
in time units. As usual G is Newton's gravitational constant 
and c is the velocity of light. The orbital inclination derived 
from ? is i = (77.4 ± 0.4)° or (102.6 ± 0.4)°. Given the 
mass function /, these values imply a pulsar mass m p — 
(1.67±O.O8)M , a total binary mass M t = (2.70 ± 0.11)M© 
and a mass ratio R s = 1.62 ± 0.03. 



3.3.2 Estimating the uncertainties 

The uncertainties quoted above were estimate d usin g the 
Bayesian technique described by ISplaver et al] l|2002l ). We 
assume that m c and cos i have constant a prion probability. 
For each point in a grid of m c , cosi values we calculate 
(trivially) the r and s parameters to describe the Shapiro 
delay, assuming that general relativity is the correct theory 
of gravitation. 

We then fit a timing solution similar to that of Table 1 
to the TOAs but keep (r, s) fixed. We record the resultant 
X 2 and calculate the probability density for that (m c ,cosi) 
space using 



p(m c , cos i) — exp 



x 2 ( 



m c , cos 1 



(2) 



where Xmin i s the minimum value of \ 2 m the whole map. 
The contour levels that include the 68.3, 95.4 and 99.7% of 
all probability are displayed in Fig. [8] as the thin contours. 
We then project the 2-D probability density function (PDF) 
into the m c and cos i axes to calculate lateralised 1-D PDFs. 

We translate the 2-D PDF of the (cosi,m c ) space to a 
2-D PDF of the (mp,m e ) space using the mass function / 
lLorimer fc Kramerl (20051 ). This is then lateralised onto the 
m p axis, resulting in a 1-D PDF for the pulsar mass. It is 
from the 1-D PDFs that we derive the uncertainties for cos i, 
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Figure 6. Top: Radio frequency of TOAs as a function of orbital 
phase. The orbital phases with the larger number and frequency 
spread of observations are those for which the DM measurements 
should be better. These are best at an orbital phase of 0.95, where 
superior conjunction happens. Bottom: DM as a function of or- 
bital phase. All TOAs were divided in 10° orbital phase bins and 
a DM was derived for them, keeping all other parameters fixed. 



3.4 Orientation of the orbit and component 

masses 

3.4-1 Kinematic effects: orbital orientation 

The additional, previously undetermined parameters in Ta- 
ble 1 are the proper motion (with total magnitude \i and 
position angle M ) and the apparent variation of the size of 
the orbit, x a ; this is given by 

io , . . n , . ti 2 D + a t x s a 

— = — ucotisin(0„ — S2) + - 1 1 — , (3) 

x c x a 

where i is the orbital inclination, Q is the position angle of 
the line of nodes, x s is the variation of x due to an intrinsic 
variation of i, a is a variation of the semi-major axis a, D is 
the distance to the pulsar and a; is the difference of Galactic 
accelerations of the binary and the solar system, projected 
along the line of sight. 




0.2 0.4 0.6 0.8 1 
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Figure 7. Top: Post-fit TOAs versus orbital phase. The Shapiro 
delay was not taken into account, but all Keplerian parameters 
were fit. Bottom: Post-fit TOAs versus orbital phase, with Shapiro 
delay taken into account. The different colours indicate different 
observing systems, as in Fig. [3] 



The first term is by far the largest; it is an effect of 
the proper motion, which constantly changes the viewing 
geometry. If this causes a change of i as viewed from Earth 
then there will be a s ecular change of the projec t ed size of 
the or bit, x = a sin i/c lArzoumanian et all (1 19961 ) ; fkopeikirj 
i| 19961 ). 

The second term describes the changing Doppler shift 
of the binary relative to the Solar System barycentre. For 
a nominal DM distance of 6.4 kpc, this amounts to —4.6 x 
10 _20 s _1 . This is three orders of magnitude smaller than the 
measurement error; this term is thus ignored in the following 
discussion. The third term is due to any real changes in the 
orbital plane of the binary system. The only likely contribu- 
tion t o this is from spin -orbit coupling ISmarr fc Blandfordl 
( 1976); ILai et al.1 (Il995l ); the resulting x s should be one or- 
der of magnitude smaller than the measurement error (see 
detailed discussion in § I3.4.4[I , Finally, the fourth term is due 
to a real change in the size of the orbit. The prediction for 
a/a caused by gravitational wave emission is of the order of 
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many orders of magnitude too small. Therefore, 



only the first term is likely to give a significant contribution 
to the observed x. 

In this case, for each of the two possible values of i 
discussed above there are two possible (Q M — fi) solutions; 
these are depicted in Fig. [5] by the intersection of the x and 
? lines. However, because of the uncertainty in the measure- 
ments of x and these two pairs of solutions merge into two 
wide areas centred at (0 M — fi) ~ 90° for i ~ 102.6° and 
(6 M -fi) ~ 270° for i ~ 77.4°. 



3.4-2 Kinematic effects: apsidal motion and masses 

The observed apsidal motion cb = (86.38 ± 0.08)"/century 
is, in principle, a combination of the relativistic apsidal ad- 
vance (jj r ; a kinematic contribution uik due to the changing 
viewing geometry of the system; and a contribution caused 
by spin-orbit coupling u) s : 



(4) 



We will now consider each of these terms in detail. 

The ojk is caused by t he change of vie wing geometry 
due to the proper motion \l iKopeikinl (|l996h : 



d>k = /U cos(6 M — fi), 



(5) 



where // = /i/sini = O'.'56/century. Given the aforemen- 
tioned uncertainty of 0,j — fi it is impossible at the moment 
to know the exact value of UJu, but it has a very sharp upper 
limit given by | cos(Q^ — fi)| < 1. This means that uncer- 
tainty in Co r caused by the lack of precise knowledge of ujk is 
at most equal to ±/i' = O'.'56/century. This is more than six 
times larger than the current uncertainty in the measure- 
ment Of LOo- 

This affects the precision of our estimate of the total 
mass of the system. If u T = Cj , then the total system mass 
could be determined directly using 



3n 



5/3 (T M t ) 2 / 3 



(6) 



where rib = 2n/Pb is the orbital frequency and Pb is the or- 
bital period (e.g. Lorimer & Kramer, 2005); the result would 
be M t = (2.697±O.OO4)M . Because of the uncertainty of to k 
and uir = L0 o — dJk, we obtain instead Mt = (2.697±0.029)A/o 
(99.7% confidence limit, see § [gX3 - 



3-4-3 Further uncertainty estimates 

To estimate these uncertainties rigorously, we have made 
a second \ 2 map that uses the same principles described 
above, but includes a third dimension, the difference in the 
position angles of the proper motion and line of nodes B,, — 
fi. This map makes use of the extra information provided 
by the measurements of x and u). 

As before, from the values of m c and cos i at each point 
we calculate r and s. But now we also use the mass function 
/ to calculate m p . From the total mass M t — m p + m c we 
calculate u r using eq.[6] From the values of cos i and (0 M — fi) 
at each point we calculate x using eq. [3] and u)k using eq. [5] 
This is made using only the best value for the proper motion 
fi, which has a relatively small uncertainty. 

We then fit the resulting timing solution to the TOAs, 



but keep the four computed parameters fixed (r, s, x and 
<jj — uj r + cjfe). We record the resultant % 2 arm calculate 
the probability density for that point of the (m c , cosi, 
O m — fi) space using an expression similar to eq. [2] This 
3-D "cube" of probability density is then projected onto its 
faces: the (cosi, m c ) space (see thick contours in Fig.[SJl and 
the (cosi, M — fi) and (m c ,0 M — fi) spaces (Fig. We 
then project this 3-D PDF onto the three axes, calculating 
lateralised 1-D PDFs for m c , cosi and O m — fi. It is from 
these that we derive the uncertainties of these parameters. 

In the (m c , O m — fi) panel of Fig. [9] we can see that if 
O fl — fi were more precisely known, we would have a much 
better estimate of m c . This is the graphical demonstration 
that the kinematic effects are now the main source of uncer- 
tainty in the measurement of m c and m p , not the uncertainty 
intrinsic to the measurement of u>. We can also see that as 
the precision of x improves, we will have a total of four pos- 
sible values for O m — fi, two for each possibility of cosi. The 
two solutions closer to the centre of Fig. [9] (0 M — fi = 180°) 
will have a large value of m c (and m p ), since, according to 
eq. [S] for this orbital orientation Wk < 0, therefore uj r > uj . 
The opposite will be true for the two solutions closer to 
O m — fi = 0° . Unless we can independently determine cos i 
or Q p — fi, we will have two degenerate values for m c and m p . 
In principle, this degeneracy can be lifted if we can measure 
the orbital annual parallax. This should be possible with a 
fourfold improvement in timing precision. 



3.4-4 Contribution from spin-orbit coupling 

Thus far we have assumed that the spin-orbit contributions 
to ui and x ar e negligible. Comparing eq. 7.54 a nd fol l owin. 
m IWilll (|l993h with the more general equation in Wex ( 199. 
we find for uj s : 



371-J2 



R c 



' 3 2 

1 sin 6 + cot i sin 9 cos 8 cos $0 



where a is the semi-major axis of the orbit of the companion 
as seen from the pulsar, R c is the companion radius (which 
we assume to be one solar radius), Pi, and e are as listed 
in Table 1, J2 is the quadrupolar moment of the star, 6 
is the angle between the orbital and companion's angular 
momentum and $0 is the longitude of the ascending node in 
a reference frame defined by the total angular momentum 
vector. For PSR J1903+0327, we obtain a maximum value 
of cus = ,h x 7.9 x 10 5 "/ century. 

As we have shown in § 13.11 the companion star has a 
mass, age and temperature similar to t hat of the Sun. Suc h 
stars should rotate slowly like the Sun llrwin et all (|2009h : 
the latter has a rotational velocity of v ro t,0 = 2kms _1 . This 
should result in a quadrupolar moment similar to that of 
the Sun, J2,q ~ 1.7 x 10~ 7 . For this value of J2 we get 
uj 3 ~ 0'.'013/century, about seven times smaller than the 
measurement uncertainty of ui B . 

However, because the companion's rotation is not well 
constrained by the optical measurements (which indicate 
Vrotsini, < 140 kms -1 , 3-a, see § 12. ip . we can use the 
agreement between /13, ? and uj r to constrain u}„ assuming 
that GR is the correct theory of gravity. The minimum total 
mass compatible (at the 1-cr level) with /13 and ? is 2.59Mq. 
The corresponding minimum ui r (84.1" /century) can then be 
used to derive a 1-er limit of uj s — 0J o —<jj r ,min < 2.3"/century 
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Figure 8. Companion mass as a function of cosi and m v . The thick contour levels are derived from a 3-D x 2 map of the cos i, m c and Q 
space (where Q is the position angle of the line of nodes) and then collapsed on the planes represented in the figure (see text for details). 
The thin contour levels represent a 2-D \ 2 map of the cosi — m c space calculated taking only the Shapiro delay into account. The lines 
represent the constraints derived from the spectroscopic mass ratio (R, black dashed), the apsidal motion (u>, solid orange - here with 
increased uncertainty due to the proper motion, solid pink), the harmonic amplitude (/13) and harmonic ratio (?) of the Shapiro delay 
(in purple) and finally an upper limit on the inclination given by x (pink solid line). The gray area in the mass-mass diagram is excluded 
by sini < 1. In the marginal plots we can see that the 1-D probability distribution functions for the pulsar and companion masses are 
much narrower when the apsidal motion (even with uncertainty caused by the proper motion) is taken into account (thick lines), but this 
assumes that there are no significant classical contributions to ui. The latter must be < 2.3"/century (1-cr) given the agreement between 
the /13, ? and uj ± fi' bands. 



(here we are assuming a median expected value of for ijjk). 
Future optical/near infrared measurements might be able to 
better constrain v Iot . If this is small then the agreement of 
the 3 different PK parameters can be used directly as a test 

of general relativity. 

To calculate x s , we use eq. 81 of I Wed l| 19981 ) to derive 

x R 2 

T = 3?F J2 P„a 2 (l-e 2 ) 2 ° 0t * C ° S ° Sin 9 Sin $ ° ' (8) 

Using the parameters above, we obtain |i s /a;| < 2.2 x 
10' 18 s"\ which is one order of magnitude below the mea- 
surement error. Therefore, an improved measurement of x 
will likely improve the determination of the orbital orienta- 
tion. 

3.4-5 Component masses and orbital inclination 

From the 3-D map, we obtain the following parameters: 
i = (77.47 ± 0.15)° or (102.53 ± 0.15)° (1-cr), M t = 
(2.697 ± 0.029) M , m c = (1.029 ± 0.008) M Q , m p = 
(1.667 ± 0.021) M© and R = 1.620 ±0.008. The underlying 
probability distribution functions for these parameters have 
a shape that is very different from a Gaussian curve (see 



Fig-©! therefore it is not very meaningful to refer to a lim- 
its. For all the previous values we indicated instead 99.7% 
confidence limits. These values and their uncertainties can 
be understood qualitatively as the result of the intersection 
of the l-<7 bands of /13 and uj t . Their intersection results in 
measurement of i that is more precise than the value derived 
from s (or Such a result cannot be understood using the 
regular r, s parameterization of the Shapiro delay. 

These values are entirely consistent with those derived 
solely from the Shapiro delay in § 13.31 but much more pre- 
cise, despite the uncertainty introduced by uik- 



4 IMPLICATIONS 
4.1 Neutron star mass 

The recent mass measurement of PSR J1614— 2230 (m p = 
(1.97 ± O.4)M , 1-cr, see Demorest et al. 2010) rules 
out the presence of hyperons, bosons and free quarks at 
densities comparable to the nuclear saturation density, 
and demonstrates that NSs can be stable at masses well 
above the Chandrasekhar mass. The mass measurement of 
PSR J1903+0327 (m p = (1.667 ± O.O21)M , see § Gy} 
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Figure 9. Confidence contours for cos i and m c as a function of the orientation of the line of nodes relative to the direction of the proper 
motion © M — f2. The contour levels include 66.3, 95.4 and 99.7% of the total probability, which is derived from a 3-D x 2 map of the 
(cosi,rfi c , M — Q) space and then collapsed onto the planes represented in the figure. The dashed contour lines indicate 90° < i < 180°. 
The vertical purple dashed lines are the constraints introduced by the harmonic ratio of the Shapiro delay (<;), the pink solid curves are 
the constraints derived from the measurement of x. In the marginal plots are displayed the 1-D probability distribution functions. We 
can see that there is still a large uncertainty in (0 M — f2). This introduces a large uncertainty in the mass estimates (of which m c is 
shown in the plot). 



supports the latter conclusion, being also inconsistent with 
the softest proposed equations of state for super-dense mat- 
ter. This plus t he ma ss measurements of PS R J0437— 4715 
IVerbiest et all (|2008l ) and PSR J0621+1002 ISplaver et all 
(2002) confirm beyond doubt the results of a previous sta- 
tistical anal ysis of masses of sev eral binary MSPs in globu- 
lar clusters iFreire et all (|2008al ). which showed that MSPs 
have a much wider mass distribution than observed among 
the members of double neutron star systems. The likely 
cause is the accretion episode that s pun up t hese NSs 
Bhattach arva fc van den Heuvell l|l99^ ; lLorirner] (|200Sh . 

There are several reasons why PSR J1903+0327 ap- 
pears to be a normal MSP (i.e., spun up by accretion of 
mass from a companion) : a) its spin period is typical of that 
of MSPs and ~8 times shorter than that of the fastest young 
pulsar, PSR J0537— 6910, near the N157B supernova rem - 
nant in the Large Magellanic Cloud iMarshall et alj (|l998l ) 
b) its magnetic field (Bo) is similar to those observed among 
MSPs and two orders of magnitude smaller than the small- 
est Bo thus far observed for a young pulsar, 3 x 10 10 G for 
PSR J1852+0040, a 105-ms X-ray pu lsar near the centre 
of the supernova remnant Kesteven 79 Hal pern fc GotthelJ 
(|201Ct) c) its mass suggests it accreted mass from a for- 
mer donor star. Precise mass measurements for compo- 
nents of double neu tron star s ystems range from 1.2 5 Mp, 
iKramer et ail (|2006T ) to 1.44 M IWeisberg etall (|2010T l; this 



shows that a large major ity of NS s form with masses near 
the Chandrasekhar limit lLorime 3 Hil). If this was also 
the case for PSR J1903+0327 then it must have gained 
AM ~ 0.22 — 0.42 Mq during accretion, which is well within 
the range expected for low-mass and intermediate-mass bi- 
nary pulsars (Pfahl et al. 2002). All of this is inconsistent 
with scenarios where PSR J1903+0327 forms directly from 
the collapse of a single star. Even alternative MSP formation 
scenarios (like merger of massive WDs or accretion induced 
collapse) require two stars to form the MSP. 

4.2 PSR J1903+0327 is not a triple system 

The discovery that the companion to PSR J1903+0327 is 
the MS star previously suspected of being associated with 
the system immediately rules out the detailed triple system 
scenario proposed in Champion et al. (2008). Our measure- 
ment of the variation of eccentricity had already ruled out 
the possibility that the Kozai mechanism is the origin of the 
observed eccentricity (Gopakumar, Bagchi & Ray 2009). 

Despite this, it is clear that the present companion can- 
not be the star that recycled the pulsar. If it were, then it 
must have at some point filled the region where matter will 
remain bound to it (its "Roche lobe") in order to transfer 
matter to the companion. The problem is that at the closest 
approach between the stars the companion is ~ 23 times 
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smaller than its present Roche lobe. Being an unevolved MS 
star it was even smaller in the pas10 and even less likely to 
have supplied the material that recycled this pulsar. 

Could the former donor have exchanged orbits with the 
main-sequence companion and now orbit the J1903+0327 
binary system at a large distance? If there were a third 
component in the system, the J 1903+0327 binary system 
would be accelerating towards it, with a line-of-sight com- 
ponent given by a;. This should produce a variation of its 
orbital period given by Pt/Pt = ai/c. Variations of a; due 
to the relative motion of the J1903+0327 binary system and 
the outer component should produce a variation of the first 
derivative of the spin frequency (i>), v (e.g., Backer, Foster 
& Sallmen, 1993). Our measurements of Pt, and v are not 
statistically significant (see Table[T|, i.e., we detect no third 
component in this system. If the former donor still exists, it 
is no longer bound to this system. 



4.3 Motion of PSR J1903+0327 in the Galaxy 

A possible explanation for the absence of the former donor 
is that it might have been exchanged by the present com- 
panion. Before the discovery of PSR J1903+0327, this 
was the only mechanism known to produce MSP bina- 
ries with eccentric orbits. This can only occur in envi- 
ronments with hig h stellar densi ties such as the cores of 
globular clusters iPhinnevI |l992) or, hypothetically, the 
Galactic Centre. This is the reason why all other ec- 
centric binary MSP s are observe d in g lobular clusters 
iRansom et all (l2004h; iFreire et al.l l|2004l ): [r 
(l2005h : I Frcirc et al. (2008b). For this reason it was proposed 



Champion et all (120081 ) that PSR J 1903+0327 might 



have formed in such an environment and subsequently have 
been ejected by the recoil produced by the presumed ex- 
change interaction. 

To investigate this possibility, we used the radial veloc- 
ity 7 and the pulsar's Galactic coordinates, DM distance, 
and proper motion to derive all three coordinates of its 
position in space and all three components of its velocity 
vector, as in Wex et al. (2000) and Lazaridis et al. (2009); 
the results are displayed in Table 1. We can thus calculate 
the past trajectory of the binary in the Galactic potential 
iKuiiken fc Gilmorel (|l989h . We implement this process in a 
Monte Carlo simulation of 10000 orbits in the model Galac- 
tic potential. We take the uncertainties of fx, 7 and distanc^B 
into account by integrating many orbits with random initial 
conditions having a distribution of parameters consistent 
with the observed values and uncertainties. We also com- 
pare our results with those obt ained usin g an a lternative 
model for the Galactic potential IPaczvriskil (|l990l ). 

Our results are presented in Figs. HOI and llll The start- 
ing orbits have a range of proper motions and systemic ve- 
locities consistent with the uncertainties in Table 1. We also 
assume a 1-cr relative uncertainty of 15% in the distance. 
For each point we integrate the pulsar orbit back in time 



for about 300 Myr and record the maximum (R max ) and 
minimum (R m i n ) distances to the centre of the Galaxy. The 
diagonal black line shown in the top plot of Fig. [TT] indi- 
cates circular orbits and the gray line eccentricities of 0.2. 
Rmax cannot be smaller than the minimum possible distance 
to the Galactic centre at present (blue line), this is given 
by 7?osinZ, where Ro is the Sun's distance to the Galac- 
tic centre (7.7 kpc in the model used for this simulation 
IKuiiken fc Gilmorel {198$ )) and I is the pulsar's Galactic 
longitude, 37.33° . In this plot we can see that the minimum 
possible distance from the Galactic centre is always larger 
than 3 kpc. This excludes formation in the Galactic centre 
and in the bulge globular clusters. 

It is also apparent from our simulations that the dis- 
tance from the Galactic plane never exceeds 270 pc (99.7% 
confidence limit), which is unusual for MSPs. The globular 
clusters observed o utside t h e bul ge have a wide distribution 
of Galactic heights [Harris! (|l996l ). which implies that most 
of them orbit the Galaxy well outside the plane. The for- 
mation of PSR J1903+0327 in such a cluster would require 
that the system is ejected at the exact time the cluster is 
crossing the Galactic plane and that the resulting velocity 
is very nearly in the Galactic plane and close to the average 
local Galactic rotation (the system's current peculiar veloc- 
ity is 37+9 km s , see bottom plot of Fig. lll[) . We find this 
possibility highly unlikely. The probability of formation in 
an exchange interaction appears therefore to be extremely 
small. 



5 FORMATION OF PSR J1903+0327 

To summarise, the measured mass of PSR J1903+0327 (plus 
its small spin period and magnetic field) make any scenarios 
where it forms as it is unlikely; all pulsar characteristics are 
consistent with it having been recycled (§ 14. ip . Our limits 
on the variation of the eccentricity and our optical detection 
exclude the triple system scenario proposed in Champion et 
al. (2008, S I4.2JI . In particular, we find that the present com- 
panion is not the donor star that recycled the pulsar; if this 
former donor still exists it is no longer bound to the system. 
Furthermore, our calculations of the 3-D systemic motion 
(§ 14.30 show that it is unlikely that the system originated 
in an environment with high stellar density. Therefore, it is 
unlikely that the former donor star was exchanged by the 
present companion. So how did PSR J1903+0327 form? And 
what happened to the former donor star? 



5.1 Multiple star spiral-in 

We suggest that the system started as a triple (indicated 
below in square brackets), where the outer companion (the 
present MS star companion) was initially on a much wider 
orbit than at present. The inner binary system (indicated 
in parentheses) had a much shorter orbital period and con- 
sisted of two somewhat more massive MS stars: 



3 If this star had ever filled its Roche lobe the orbit would have 
been circularised. This is inconsistent with the present large and 
non- varying eccentricity. 

4 We used a 1-tr relative uncertainty in the distance of 15% such 
that it is similar to the range of optically derived distances. 



[(MS + MS; 12 + 2M Q ;2Od)+MS;lM ;lOOOd] 

(the numbers are conjectural, we show them for illustrative 
purposes only). The more massive star (the progenitor of 
the pulsar), evolves and becomes a red supergiant (RSG): 
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Figure 10. Present Galactic positions (dots) and past orbits, integrated for 1 Gyr (solid curve s) of the four MSPs with known 
radial velocities: PSR J1 9 03+0 3 27 (Red), PSR J1012 +5307 [Dark Blue, lLazaridis etafl (2009h ]. PSR J1023+0038 [Light Blue, 
iThorstensen fc Armstrong! J2005I) ; lATchibald et all i2009h l and PSR J1738+0333 [Olive Green]. The Sun's position is indicated by the 
yellow circle. In the top plot we look at the Galaxy from the North Galactic pole, in the bottom plot from the direction in the Galactic 
plane perpendicular to the line from the Sun to the Galactic centre. We can see that PSR J1903+0327 is always near the Galactic plane, 
and that its orbital velocity about the Galaxy has a relatively small eccentricity. 



[(RSG + MS; 12 + 2M ; 20d) + MS; 1M ; lOOOd] 

At this point, it starts transferring mass to the other compo- 
nent in the inner binary. In such a situation, mass transfer 
is generally unstable, and the envelope of the RSG is ex- 
pected to engulf the companion and lead to a common en- 
velope (CE) phase, where the helium core of the RSG and 
the inner companion are embedded in the envelope of the 
RSG. Because of friction with the envelope, the orbit of this 



embedded binary will decay, releasing orbital energy in the 
process. If this energy is sufficient to eject the envelope, this 
will leave a much closer binary consisting of a helium star 
and the more-or-less unaffected inner companion (see, e.g., 
Bhattach arva fc van den Heuvellll99ll ): 

[(He + MS; 3 + 2M ; 3d) + MS; 1M Q ; lOOOd] 

If the outer component of the triple system is close enough, 
it may also be affected by this CE phase, since the en- 
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Figure 11. Monte-Carlo simulation of 10000 pulsar orbits. Top: 
Minimum and maximum distances from the centre of the Galaxy. 
We can see here that the pulsar is never within 3 kpc from the 
centre of the Galaxy Bottom: Starting velocities in previous simu- 
lation relative to each simulated pulsar's average of local velocities 
(i.e., the system's peculiar velocity). 



velope will greatly expand because of the orbital energy 
that is b eing deposited within it (to radii of ~ 5 — 10 AU; 
see, e.g.. |Podsiadlowskll200ll ). It will also be engulfed by 
the e xpanding envelope and will also experience a spiral-in 
phase lEggleton fc Verbunt] (1 19861 ). This will produce a much 
closer triple system. Without realistic hydrodynamical sim- 
ulations, it is difficult to determine the post-CE parameters. 
Here, we will consider two potential outcomes, which we will 
discuss in more detail later: 

a) [(He + MS; 3 + 2M ; 2d) + MS; 1M ; lOOd] 

or 



Eventually, the He core will explode in a supernova and pro- 
duce a NS. The sudden mass loss associated with this event 
will make the orbits of the other two components wider and 
somewhat eccentric: 

a) [(NS + MS; 1.4 + 2M ; 2.5d) + MS; 1M©; 120d, e = 0.1] 
or 

b) [(NS + MS; 1.4 + 2M ; 8hr) + MS; 1M ; 70d, e = 0.44] 

This event could be either a normal Fe-core collapse super- 
nova (SN) or an electron-c apture (e-capture) supernova (see 
iPodsiadlowski et al.l 120(53 . and further references therein). 
Given the fact that the system remained bound after the ex- 
plosion, it is not likely that the SN produced a major kick, 
nor large fractional mass loss, particularly under scenario 
"a", where the system is wider and acquires a small eccen- 
tricity. This is consistent with the observed peculiar velocity 
for the system, which is smaller than for the other systems 
where we can make a complete determination of this quan- 
tity (see Fig. I10p . We note too that such SN characteristics 
are expected from e-capture supernovae. 

5.2 Post-supernova evolution 

At this stage the inner binary resembles the progenitor of a 
typical low- or intermediate- mass X-ray binary (L/IMXB). 
Once the normal stellar component of the inner binary has 
evolved to fill its Roche lobe, mass transfer to the NS will 
start to increase its mass and spin it up to millisecond pe- 
riods (see, e.g., Pfahl et al. 2002). Initially, when the sepa- 
ration of the inner binary is sufficiently small compared to 
the periastron of the outer component, the evolution of the 
L/IMXB will not be significantly affected by the presence of 
the outer component. 

If the mass-transfer rate is very high, mass loss from the 
inner binary will cause the orbit of the outer component to 
widen to some degree. However, even for conservative mass 
transfer the inner binary is generally expected to widen more 
drastically, in particular in the 2.5 day scenario (scenario "a" 
in § 15.11 see Podsiadlowski et al. 2002): 

[(MSP + MS; 1.67 + O.3M ; 50d) + MS; 1M ; 150d, e = 0.1] 

until its separation reaches a critical fraction of the perias- 



b) [(He + MS; 3 + 2M ; 0.8d) + MS; 1M ; 50d] 



PSR Jl 903+ 0321 15 



tron distance of the outer compan ion. This leads to a chaoti c 
three-body interaction (see, e.g., lHulll984l ; IPhillipsI Il993h . 
One of the most probable outcomes is the ejection of the 
least massive component - in this case most likely the mass 
donor in the inner binary, leaving a bound, eccentric binary 
system consisting of the MSP and the formerly outer com- 
panion, in an orbit tighter than before: 

(MSP + MS; 1.67 + 1M ; 95d, e = 0.44) 

If the starting orbital period is small (Pb < 8 hr, scenario 
b in § 15. then gravitational radiation might lead to the 
destruction of the remnant of the former donor by the newly 
formed MSP, very much in the same way as PSR 1957+20, 
the black- widow pulsar (Fruchter et al. 1988), is evaporat- 
ing its companion (although in the latter case the process 
might take several Hubble times). This is a possible forma- 
tion channel for isolated MSPs. Since these represent ~20% 
of the known MSPs in the Galactic disk, this should not be 
too unlikely a scenario. The end result would also resemble 
the J1903+0327 binary system . 

5.3 Alternative scenarios 

One can imagine alternative triple scenarios to the one out- 
lined above. For example, it is possible that the progenitor 
of the MSP was initially the outer component of a stable 
triple system. When this massive star evolved to become 
a red supergiant, it could have engulfed the inner binary 
completely, leading to a spiral-in of that binary inside the 
common envelope. As it spirals in, its orbit could become 
disrupted by the tide induced by the more massive compo- 
nent; both of its components then start orbiting the core 
of the more massive He star independently, but at differ- 
ent distances given their different orbital velocities prior to 
disruption. This would produce a result quite similar to the 
double star spiral-in. Another possibility is that this binary 
survived the CE phase, but it was later disrupted when its 
orbital radius increases due to a mass transfer stage. 



6 CONCLUSIONS 

Our Arecibo and Green Bank radio timing observations in- 
clude a full determination of the relativistic Shapiro delay, 
a very precise measurement of the apsidal motion and new 
constraints of the orbital orientation of the system. Through 
a detailed analysis of all of these, we derive new constraints 
on the masses of the pulsar and companion. The mass of 
the pulsar (1.667 ±0.021 Mq) adds to a growing population 
of NSs with masses which significantly exceed 1.4 Mq and 
rules out some of the proposed equations of state of super- 
dense matter. This high mass also suggests that the pulsar 
was recycled. 

Our observations with the Very Large Telescope reveal 
shifts in the spectral lines of the suspected companion to 
PSR J1903+0327, which show that it has the 95-day orbital 
motion expected for the pulsar's binary companion. This 
star has characteristics very similar to the Sun; making this 
system unique among all known binary MSPs. A detailed 
discussion of the results of the radio timing and optical spec- 
troscopy rules out most proposed formation mechanisms. 



We are then forced to accept the only remaining possibility, 
i.e., that the system originated from a triple system with a 
compact inner component and the presently observed MS 
companion as the outer component. The inner compact bi- 
nary then formed the observed MSP, with the elimination of 
the former donor star (i.e., evolution towards an "isolated" 
MSP), or its ejection from the system. 

Detailed stellar evolution simulations of the formation 
of this system might provide new insights on how MSPs 
form. As an example, the mass loss and kick velocities asso- 
ciated with the supernova event that formed this NS must 
have been small enough to allow the whole system to remain 
bound. Such a gentle formation would produce a small re- 
coil velocity for the whole system; this is consistent with the 
small peculiar velocity observed at present. This would gen- 
erally favour alternative small-kick formation mechanisms, 
like electron capture supernovae. This means that a dynam- 
ical study of this system could lead to a conclusive demon- 
stration of alternative formation channels for NSs in general. 
Some of these exotic mechanisms (for example, white dwarf 
mergers) could be especially important for explaining the 
existence of isolated MSPs. 
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